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Abstract 

The spectral properties and phase diagram of the exact integrable spin one quantum chain intro- 
duced by Alcaraz and Bariev are presented. The model has a U(l) symmetry and its integrability 
is associated to an unknown R-matrix whose dependence on the spectral parameters is not of dif- 
ference form. The associated Bethe ansatz equations, that fix the eigenspectra, are distinct from 
those associated to other known integrable spin models. The model has a free parameter t p . We 
show that at the special point t p = 1 the model gets an extra £7(1) symmetry and reduces to the 
deformed SU(3) Perk-Schultz model at a special value of its anisotropy q = exp(i2ir/3) and in 
the presence of an external magnetic field. Our analysis is done either by solving the associated 
Bethe-ansatz equations or by direct diagonalization of the quantum Hamiltonian for small lattice 
sizes. The phase diagram is calculated by exploring the consequences of conformal invariance on the 
finite-size corrections of the Hamiltonian eigenspectrum. The model exhibits a critical phase ruled 
by c = 1 conformal field theory separated from a massive phase by first-order phase transitions. 
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I. INTRODUCTION 



The anisotropic spin 1/2 Heisenberg model, or XXZ quantum chain, and the 6- vertex 

n 

model are considered as paradigm of exact integrability in statistical mechanics [lj. In the 
XXZ quantum chain the ^-component of the total magnetization is a good quantum number 
(£7(1) symmetry). Its simplest integrable generalizations that keeps the U(l) symmetry are 
spin-1 quantum chains. Models on this class are the Fateev-Zamolodchikov model [2J, the 
Izergin-Korepin model [3|, the supersymmetric OSP(l/2) model 4] and the biquadratic 
model (5)]. The integrability of these models is a consequence of the existence of a known 
associated R-matrix satisfying the Yang-Baxter equation. The associated R-matrix for these 
models are regular, i.e., depends only on the difference of the spectral parameters. 

A new exact integrable spin-1 quantum chain was derived by using the coordinate Bethe 
ansatz jsj, or a matrix product ansatz Q. The derivation of the integrable model through 
these last approaches does not depend on the knowledge of the associated R matrix. Dis- 
tinct from the other integrable and non integrable [8| spin-1 models, whose physical 
properties are well studied, almost no physical information is known for this new quantum 
chain, besides its exact integrability. The unknown associated R-matrix is not regular |g| 
since it does not satisfy the Reshetikihin criterion [9j. The Bethe ansatz equations (BAE) 
that fix the eigenenergies are also quite distinct from the corresponding equations of other 
spin-1 integrable quantum chains. 

In this paper we are going to present an extensive analytical and numerical analysis of 
the eigenspectra properties of this new spin-1 quantum chain. Based on solutions of the 
associated BAE, whenever it is possible, and diagonalizations of the quantum Hamiltonian 
on small lattices (L = 2 — 24) we are able to predict some of its critical properties. 

The paper is organized as follows. In section 2 we present the model and the BAE that 
fix the eigenspectra. The model has a U(l) symmetry and is exact integrable for any value 
of a free parameter t p . We show, that for the special value t p = 1, the model is related to 
the deformed SU(3) Perk-Schultz model with deformation parameter q = exp(i27r/3) lo| 
in the presence of an external magnetic field. In section 3 we analyze the eigenspectra of 
the quantum chain in several regions with distinct values of the free parameter t p . Based on 
conformal invariance predictions the critical properties of the model are calculated. Finally 
in section 4 we summarize our results and present our conclusions. 
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II. THE MODEL 



Instead of presenting the quantum Hamiltonian in terms of spin-1 SU(2) matrices 
(S x , S y , S z ) it is more convenient to present it in terms of the 3x3 Weyl matrices E l,m 
(l,m = 0,1,2), with elements E l,m )ij = 5i t i5 m j. At each lattice site i we may have zero 
particle (n« = 0), one particle (rii = 1) or two particles [rii = 2) or equivalently Sf = — 1, 
Sf = and S* — 1, respectively. The dynamics of these particles, in a periodic chain with 
L sites, is described by the Hamiltonian 



where 



J = l a,^,7,/i=0 j=la=l 



pl,0 _ F 0,l _ p2,l _ pl,2 _ 1 F 0,2 _ F 2,0 _ , 
1 0,1 — 1 1,0 — 1 1,2 — 1 2,1 — — 1 > 1 2,0 — 1 0,2 — ~ l Vi 



r ; 2 = = e-^v^T, r$ = r? : ° = -e^^T, (2) 

are the hopping parameters (off diagonal) and the static terms (diagonal) are given by 



F o,o _ F 2,2 _ n F i,o _ r o,i _ L 
1 0,0 _ 1 2,2 — u ; 1 1,0 _ 1 0,1 _ ^ ' 

r i.i _ , ' F o,2 _ F 2,o _ tp 

1 1,1 — l p 1 2^~> 0,2 — 1 2,0 _ 2" , 

pl.2 _ 1 , .V3, p2,l _ 1 .V3, /o\ 

1 1,2 — + i ~^ l py 1 2,1 — ^jj" i ~^ l p- \ 6 ) 

The parameter i p is free and /i plays the role of a magnetic field in the z direction (z-magnetic 

field) or a chemical potential controlling the magnetization or the number of particles in the 

ground state, respectively. 

The Hamiltonian ([[]) is non hermitian. It is interesting to observe that its non hermiticity 

is not only due to the presence of complex matrix elements in the diagonal (see (JH])) (as 

happens in the quantum deformed SU q (N) models) but also due to the presence of complex 

non diagonal elements (see ©). The Hamiltonian ([[]) has a U(l) symmetry due to its 

commutation with the total density p of particles: 

n L 2 

p = —, n = Y, rij, rij = + 1 = ^ &Ef a . (4) 

L j=l a=l 

As a consequence its associated eigenvector space can be separated into disjoint sectors 
labeled by the the total number of particles n (or density p) or equivalently by its magneti- 
zation. 
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The quantum chain ([I]) corresponds to one of the exactly integrable models introduced 
in {(J. It is given by the choice e = 1 in Eq. (15) of 6]. As compared with the original 
presentation of the model, we also added a harmless ^-magnetic field so that the ground state 
of ([I]) at h = has, for any value of t p total density p = 1 or equivalently zero magnetization. 
The Hamiltonian ([T]) is exact integrable for arbitrary values of the parameter t p and magnetic 
field h. 

At t p = 1 the non-diagonal couplings T^'q = Tx'i = ^I'l = ^11 = and the model has an 
additional £7(1) symmetry. The number of sites with single and double occupancy are now 
conserved separately. The parameter t p can be interpreted as an anisotropy parameter, being 
t p = 1 the isotropic point. At this isotropic point, apart from a contribution Y.i=i(E^'° ~ 
that vanishes in the periodic chain, the Hamiltonian, with h = 0, is given by 

= l,h = 0) = H PS {if) - \ £ Ef + |, (5) 

where 

Hps{i) = EE {cosh -yETEgl 

j=l a =0 

+ £ [^^E^E^-E^E^ + E^E^ + ^E^m (6) 

/3=a+l j— | J 

is the deformed spin-1 SU(3) Perk-Schultz model [10| at the special value of the deformation 
parameter q = e 7 , 7 = z27r/3. This model is also known as the anisotropic 577(3) Sutherland 



model It is important to stress that the related Perk-Schultz Hamiltonian is the ferro- 
magnetic one (signal +) in the presence of a special magnetic field (value h 1 = 3/2, h 2 = 0) 
favoring single occupied sites. A simple calculation show us that the ground state for the 
related Perk-Schultz model happens in the sector with total density p — 1. It corresponds 
to the trivial state |11 ■ ■ -1 > where all the sites are single occupied. The model (JT]), for 
t p 7^ 1, can then be considered as the spin-1 anisotropic Perk-Schultz model at q = e t2w / 3 
with an additional parameter that breaks partially its symmetry. For arbitrary values of t p 
the eigenenergies and momentum are given by (| 

E = -n(— + h) + 2 cos k v P = 5>i> ( 7 ) 
zt p j=i 3=1 

where {kj = k(Xj);j = 1, . . . , n} are the roots obtained from the BAE 

e , tjL = _^ 3 mh(A,-A,-^/3) 

f-J^ smh(Aj — Xi + i2ii/ 3) 



n Si "ha~^/sv J ' do) 

" smh(Aj — A ; + «27r/3) 



with 

sinh Aj - iJ3% + (4*2 - l)sinh 2 A 7 - 

= ? — -. (9) 

t p (sinh Xj + 3 cosh Xj) 

As we see from (JH]) the left-hand side of the BAE (JSJ) are quite distinct from the corresponding 

equations for other exact integrable quantum chains. In order to consider the bulk limit 

(L — > oo) it is then necessary to study these equations for small lattice sizes. These studies, 

as we are going to see in the next section, will give us educated guesses for the topology of 

the roots {Aj} related to the low lying eigenvalues. 

Before closing this section it is interesting to consider the BAE (IE])-© at t p = 1. In this 

case they are given by 

sinh(Aj — in/3) L 
sinh(Aj + in/3) 

These equations coincide with the BAE of the XXZ chain [1] at the special value of its 
anisotropy A = — (q + 1/q) = 1/2 (q = e t2n ^ 3 ). However, in the XXZ chain the density 
of particles is restricted to < p < 1 while in the model (JTJ) < p < 2. Although the 
completeness of the Bethe ansatz solutions is always a complicated problem, we expect that 
at t p = 1 the solution obtained from ( fTUl) is not complete. Due to the additional symmetry 
(conservation of the number of pairs of particles) we should start again the Bethe ansatz jf| or 
the matrix product ansatz [7] taking into account this new symmetry. In this case we obtain 
the nested BAE of the deformed SU(3) Perk-Schultz model with the deformation parameter 
value q = e ,27r / 3 . At this point this last model is special. In 12[ several conjectures about 
the eigenspectra of this model on its antiferromagnetic regime was done. As is well known 
the eigenspectra of the anisotropic SU(3) Perk-Schultz model contains all the eigenvalues of 
the XXZ quantum chain with the same anisotropy. The equivalence ([5]) then indicates that 
the whole eigenspectra of the XXZ with anisotropy A = 1/2 is contained in the eigenspectra 
of the Hamiltonian (JTJ at £ p = 1. Moreover, a direct diagonalization of (JTJ, with t p = 1 
and h = 0, shows that for small even lattice sizes (L < 12) the low lying eigenvectors, for 
densities p < 1, coincide with those of the XXZ chain at the anisotropy A = 1/2. This 
means that there is no double occupancy of particles on these eigenstates. This can be 
understood from (J5j) due to the presence of the magnetic field favoring single occupied sites. 
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III. EIGENSPECTRA CALCULATIONS 



In this paper we restrict ourselves to the cases where the parameter t p is real and positive. 
The eigenvalues of H(t p ,h) are the same as those of —H(—t p ,—h). The Hamiltonian ([I]) 
although having a real trace is non hermitian. The exact eigenspectra calculations for lattice 
sizes L < 12 show that most of the eigenenergies are real. All the low lying eigenvalues being 
real numbers. Imaginary eigenvalues, appearing in complex-conjugates pairs, happen only 
for high excited states in the eigenspectrum. We also verify that when h = the ground 
state of ([1]) belongs, for general values of t p , to the sector with density p — 1. Moreover, for 
h = 0, the low lying excited states in the sectors with densities p = 1 + ^ and p = 1 — ^ 
(m = 1, 2, . . .) are degenerated. This degeneracy is not valid for higher excited states since 
the model (Q]) at h — does not have the symmetry under the conjugation of particles: 
<-> 2, 1 1. 

At t p = 1, where the model ([1]) recovers the deformed SU q (3) Perk-Schultz model at 
q = e l2n / 3 (see ([6])) the model is massive (non zero gap). Our numerical results indicate 
the same massive behaviour for all values of t p > as long as h — 0. As h decreases 
(negative values) it reaches the critical field h c (t p ). For h < h c {t p ) the ground state changes 
continuously its density of particles p = p{h) < 1. We do expect, in the plane (t p , p) a phase 
diagram with massless (critical) behaviour for p < 1. In this critical regime the long-distance 
fluctuations should be ruled by a underlying conformal invariant field theory. The conformal 
central charge c of the continuum theory can be estimated from the finite-size corrections of 

n 

the ground-state energy E (L,t p ,p) of the finite-size L chain [13], i. e., 

+ o(L- 2 ), (11) 



E (L,t p ,p) _ TXV S C >, 



L °" I? 

where eoo is the bulk limit of the ground-state energy per site and v s = v s (t p , p) is the sound 
velocity that can be inferred from the energy-momentum dispersion relations. Moreover, 
for each primary operator $a,A; with dimension x$ = A + A and spin sq> = A — A in 
the operator algebra of the underlying conformal theory, there exists an infinite tower of 
eigenstates, whose energies m , (L) and momentum P* m , behave asymptotically as [bJ ] 



27T 

E^ m ,(L) = E (L) + —v s (x$ + m + m') + o(L- 1 ), 
2n 

p m, m > = -£-(s 9 + m-m'), (12) 
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with m, m' = 0,1,2,.. .. 

Our spectral analysis of the Hamiltonian ([1]) was done by solving numerically the BAE 
(JHJ)-© using a Newton type method, whenever it was possible. Since there is no numerical 
method that warranties the solution for the non linear equation (jHJ), the success in finding the 
solutions depends very much on the ability to provide educated guesses for the approximated 
values. In the cases where we were not able to obtain the solutions of the BAE our analysis 
were based on direct numerical diagonalizations, or by using approximated methods like the 
power method. On these last cases our analysis were limited for lattice sizes up to L = 24 
sites. 

According to the different topologies of roots of the BAE ©-0 we divide the plane 
(t p , p) in five regions (see Fig. [Q. 




_j i i i i i i i i 

0.5 1 1.5 2 2.5 

t 

P 



FIG. 1: Phase diagram of the quantum Hamiltonian (pQ)-©. In regions 1 and 2 the ground state 
is described by real roots of the BAE. In the other regions the ground state, besides real roots, 
also contains complex ones. Regions 1,2,3 and 5 are critical and governed by c = 1 Coulomb gas 
conformal field theory. Region 4 is characterized by several crossing of the eigenenergies producing 
an oscillatory behaviour in the finite energy gaps and momentum. The lines separating regions 1 
and 5, and regions 2 and 3 are obtained by solving (|19j) - f)21[ ) and ()23p . An schematic line (heavy 
line) where we expect a first-order phase transition, separating region 3 and 4, is also shown. 

In general the roots {Xj} of (jHJ) are complex numbers. Numerical analysis, based on direct 
solutions of (jHJ), as compared with brute force diagonalizations of the quantum chains shows 
that in regions 1 and 2 (see FigjTJ) the roots {Xj} corresponding to the low lying eigenvalues 
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are all real numbers. For these cases, the right-hand side of (jHJ) is unimodular. The roots are 
then constrained or not on an finite real interval, depending on the values of t p . If t p < 1/2 
(region 1) the real roots are constrained on the interval 

- n(t p ) < Xj < n(tp), (13) 

where 

1 2t 2 + 1 

fi(g = -cosh- 1 ( r ^). (14) 

For t p > 1/2 (region 2) these roots are unconstrained. For these real roots the BAE (JED-© 
take the simple from 



Lk(X j )=2nQf + ~£<t>(\j,-\ t ), j = l,...,n, (15) 



k=l 



where 



jfe(A) = 2 (A) - Qx (A), 0(A) = 2arctan(^^; 

v3 



, /3*2 + (4t| - 1) sinh 2 A 

02(A) = — arctan(- — — ), 

smh A 

6*i (A) = arctan(\/3 coth A), (16) 

and Qf (j = 1, . . . , n) are integers or odd-half integers, depending on the particular eigen- 
state. As before the eigenenergies and momenta are given by (TH). 

In order to illustrate we present in table [TJ some eigenenergies obtained by solving the 
BAE (fl~5|) in region 1. We take in ([TJ) the magnetic field h = 0. They are obtained for several 
lattice sizes L at t p = 0.2, 0.3 and 0.4 and densities p = 0.6, 0.65, and 0.7, respectively. 

They are zero (mod. 7r) momentum eigenstates with roots {A^} symmetrically distributed 
around the origin. The corresponding quantum numbers in ( Fl6|) are Qj = ±(— ^ + j) 
(j = 1, . . . , n — pL). The finite-size corrections of the ground state energies indicate us that 
in region 1 and 2 the Hamiltonian ([Q) is critical and conformally invariant. The relation (II ip 
gives us an estimate for the conformal anomaly c of the underlying conformal field theory. 
In Fig. |2] we show, as an example, the ground-state energy per site as a function of 1/L 2 , 
for the quantum chain with t p = 0.2, density p = 0.6 and h = 0. We clearly see a linear 
behaviour as predicted by conformal invariance (see (TTTTl ) . The linear coefficient of the dot 
line in Fig. |2] gives us, from ( ITT]) , an estimate for the product 2nv s c. The sound velocity 
v s is more difficult to estimate from numerical solutions of ( lT5|) -( fT6|) . since it demands the 
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L \ (t p , p) 


(0.2,0.6) 


(0.3,0.65) 


(0.4,0.7) 


10 


-2.24745586 




-1.70579604 


100 


-2.23726600 


-1.86555023 


-1.69990939 


200 


-2.23719021 


-1.86548482 


-1.69653825 


1000 


-2.23716569 


-1.86547050 


-1.69651539 


Ext. 


-2.23716439 


-1.86546561 


-1.69624038 



TABLE I: Lowest eigenenergies per site of the Hamiltonian (JTJ) with h = for some values of L, 
t p and density of particles p. These points belong to region 1 (see Fig. [[]). They are obtained by 
solving directly the BAE (|15p . They last line is the asymptotic value obtained from the solution 
of dT9D-d22D. 

-0.7371 
-0.7372 
-0.7373 
-0.7374 
-0.7375 
-0.7376 

0.0001 0.0002 0.0003 0.0004 

L- 2 

FIG. 2: Ground state energy per site of the Hamiltonian ([IJ) as a function of 1/L 2 , for t p = 0.2, 
p = 0.6. We set h = in the figure. This point belongs to region 1 in Fig. [TJ 

calculation of excited states with non zero momenta. The BAE roots on this case are not 
symmetric. A possible way to circumvent this problem is to seek for excited zero momentum 
states related to higher values of m = m! in the conformal tower ( |12p of the identity operator 
(xcj, = 0). Our numerical analysis indicate that these energies have a configuration of BAE 
roots where pL — 2 of them are real and symmetrically distributed around the origin, and 
the two out-most roots are in the form A = ±(a + in), with a G 9ft. There is a difficulty in 
using this set of eigenlevels. As we change the lattice size their relative positions m = m! in 
the conformal tower are not fixed. However the estimates for the sound velocity obtained 
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by direct diagonalization (L < 24) of the quantum chain, although with low precision, are 
enough to indicate the position m — m' of the eigenlevel in the conformal tower. Using 
this procedure we computed the conformal anomaly at several points in regions 1 and 2 (see 
Fig. d]) obtaining c = 1.00(1). 

Since in regions 1 and 2 the ground state is described by real roots of the BAE, it is not 
difficult, in this case, to consider the bulk limit (L — > oo) of the BAE (IT5T) - (TX6|) . Following 



Baxter 



15] we define the variables 



9i 

L 



1 

2^ 



L 1=1 



(17) 



When L — )• oo, — )• s become a continuous variable in the interval x m \ n < x < x max 
satisfying the integral equation 



L 



1 

2^ 



k{X 3 



0(A(x) - X{y))dy 



(18) 



Since for the ground state the {Qj} are equally spaced by the unity, o"l(A) = will give 
us, for L — > oo, the local density of roots doc (A), that satisfies the integral equation 



where from fTTol) 



Coo (A) 



dk(X) 



1 f dk(X) 
2^ \ dX 



x/3 



2-v / 3a 00 (A / ) 



A 2cosh(2(A-A / )) + l 



dX' 



(19) 



dX 2cosh(2A) + l 



1 + 



-\/3 cosh A 



(20) 



l{Ul -l)sinh 2 A + 3^ y 

and A = X (t p , p) gives the extreme values of the roots. The total density of particles and 
the ground-state energy per site are given by 

rA 



P 



An 



cr OQ (X)dX, 



(21) 



e (t p ,p,h) = e (t p ,p,0)-hp=-p(^- + h) + 2 [ ° cos((A;(A))(T 00 (A)dA, (22) 

where h = h(t p , p) is the magnetic field that fixes the ground-state energy at the density p. 

For t p < 1/2 (region 1), X (t p , p) is always finite, i. e., X (t p ,p,h) < Q(t p ) (see f[T4"|) ). In 
Fig. El we show the ground-state energy e (t p ,p,h) for several values of t p in region 1 (we 
set h = in the figure). They are obtained by using in (j22|) the o"oo(A) obtained from the 
numerical solution of the coupled integral equations f|T9|) -f l2T|) . As we see from this figure 
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-0.72 






1 1 ' 


0.49 r 








- 0.20 








-0.74 






/ - 








^ 0.25 








-0.76 








/ 








0.30 






p 


-0.78 




0.35 


0.45 / 




Jo. 


-0.8 




0.40 




- 


U 












-0.82 












-0.84 












-0.86 


, 





i 1 i 1 






0.5 0.55 


0.6 0.65 0.7 0.75 0.8 


0.85 0.9 


0.95 



FIG. 3: Ground state energy per site of the Hamiltonian JT| as a function of the density p, for 
several values of t p in region 1. In the figure we set h = and the values of t p are shown. 
These curves are obtained by solving (|19|) - ()21|) . For a given value of t p the curves are shown for 
P < Pmax.{tp)-, with p max (i p ) given by ([23]) . For p > p max {t p ) there exist complex roots in the BAE 
solutions for the ground state and the model is in region 5. 

these equations give us a maximum total density of particles p = p mSiX (t p ) < 1 for the 
quantum chain. At the endpoints of the curves we have Ao = Q{t p ) (see f[T4|) ). Above this 
density, which we refer as region 5, some of the BAE roots associated to the ground state 
of the quantum chain have complex values. In Fig. [1] the line separating regions 1 and 5 
gives, for a given t p , the maximum density obtained from ( fT9l) - (l2Tl) . This line was obtained 
by solving numerically (jT9]) - (l2~l"j) . In order to compare with the finite-size results we give in 
the last line of table [J the estimated value of €o(t p , p, h = 0), obtained by solving ffl"9l - fl2~T]) 
for t p = 0.2, 0.3 and 0.4 with p = 0.6, 0.65 and 0.7, respectively. 

For t p > 1/2 (region 2) the roots {Xj} are not constrained and Ao = Xo(t p , n) is arbitrary. 
The maximum density compatible with only real roots of the BAE for the ground state 
is obtained by setting Ao — > oo in ( |T9l) and ( 12~TT) . In this case we can solve (|T9l)-(l211 by 
using Fourier transforms. After some long, but straightforward calculation we obtain the 
maximum density 

3^ f V3X 



Pmax 



1 + 



dX. 



(23) 



/o 2 cosh(2A) + 1 ^ ' ^(4^ _ !) sinh 2 A + 3t s 

At special values of t p we are able to solve (|2"3"|) analytically: p max (l/2) = 1, p m ax(l) = 1/2, 
Pmax(oo) = 1/4. In Fig. [1] the curve separating region 2 from 3 was obtained from the 
numerical evaluation of fl23l . As happened in region 5, for p > p ma , x (t p ) (regions 3 and 4), 
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the BAE solutions giving the ground state contain complex roots besides the real ones. 

As in region 1 (see table HJ we also calculated the finite-size corrections for the ground- 
state energy in several points of region 2. Using (llip and the same procedure as before 
our results indicate that, like region 1, region 2 is also critical and conformal invariant with 
c = 1. We then have in both regions (1 and 2) a massless behaviour with ground-state 
energy given by real roots of the BAE. This critical behaviour is quite similar as that of the 



XXZ quantum chain in the presence of a magnetic field [161. 1171]. We then expect in regions 
1 and 2 a physical behaviour described by an underlying c = 1 Coulomb gas conformal field 
theory. The anomalous dimensions of operators being given by 

xi, m = l 2 x p + m 2 /4:Xp, (24) 

with I, m = 0, ±1, ±2, . . .. The dimensions l 2 x p are obtained from (fT2|) . by considering the 
difference between the ground-state energies in the sectors with densities p and p + l/L. The 
dimensions x , m are calculated by using in ([12]) the mass gaps associated to eigenstates with 
the same density of particles. 

The parameter x p = x p (t p , p) in regions 1 and 2, which are estimated from the finite-size 



corrections of the energy (|22j) . can be calcu 
our relations ffT8|) - fT22l the method used in 
field. We obtain 



ated analytically. This is done by applying to 



161 ] for the XXZ quantum chain in a magnetic 



x p = (2£(A ))- 2 , (25) 



where £(Aq) is the dressed charge [l8| evaluated at the Fermi surface Aq of the effective 



Coulomb gas. This function satisfy the integral equation 

m = 1+ £ r- jw dA , (26) 

^ ; 7T J-x 2cosh(A - A') + 1 

In Fig. H]we show, for several values of t p in regions 1 and 2, the dimensions x p = x p (t p , p) 
obtained by using in (125]) the numerical solutions of the coupled integral equations (|T9|) . 
( 12~T|) and (1261) . We can see from this figure that for any t p , as p — > 0, x p — > 1/4. This can 
be understood from the fact that at this limit the interacting potential energy is negligible 
when compared with the kinetic energy (hopping terms). We have essentially non-interacting 
particles, where x p has the value 1/4. We also see from this figure that, for t p < 1/2 (region 
1), the limiting value of x p depends on the value of t p . This value is obtained by choosing 
A = Q(t p ) in ([26]) . where fl(t p ) is given by (|T41) . In Fig. [5] we show the limiting values of 
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Xp as a function of t p , for t p < 1/2. It varies from x p = 1/4 to x p = 1/6 as t p goes from to 
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FIG. 4: The values of x p (t p , p) (see ([24]) ) as a function of p for several values of t p in regions 1 and 
2. The curves with t p < 1/2 and t p > 1/2 belong to regions 1 and 2 of Fig. [1] respectively. The 
endpoints of the curves are the densities separating regions 1 and 2 from 3 and 5, respectively. For 
t P > 1/2 the endpoints of the curves is x p = 1/6 ~ 0.16666. 




t 

p 

FIG. 5: The values of x p at the line separating regions 1 and 5 (see Fig. [1]). The exponents varies 
from Xp = l/4tox p = l/6~ 0.1666, as t p goes from to 1/2. 

1/2. Moreover, Fig. H]shows that the limiting value is x p = 1/6 for any t p > 1/2. This should 
be the case since for any t p > 1/2 the maximum value of Ao is infinity and consequently 
( |25l) and ( 1261) give us the same result for any value of t p . The exact value x p = 1/6 can be 
understood from the relation (see Sec. II) of the model at t p — 1 and the XXZ chain with 
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anisotropy A = 1/2. When Ao — > oo the density p — > 1/2. At this density the XXZ has no 
magnetic field and its exponent is given by x p = (n — cos _1 (— A))/27r = 1/6. These results 
imply that, at the line separating regions 1 and 5, x p varies continuously from 1/4 to 1/6 
and, in the line separating regions 2 and 3, it remains fixed to the value 1/6. 

In regions 3 and 5 of Fig. [TJ some of the roots of the BAE corresponding to the ground 
state have complex values. In fact in some points of these regions we were able, for small 
lattice sizes, to solve the BAE for the ground-state energy. We verified that besides real 
roots we also have pairs of 2-strings (pair of roots of type \± = a + ib, with a,6 G SR). The 
mixture of real and complex roots in the BAE (jSJ) produces numerical instabilities turning 
the numerical solution of the BAE a quite difficult task. Due to this difficulty, instead of 
solving the BAE (JHJ) we have used the power method to calculate the lower part of the 
eigenspectrum. On this case, due to computer limitations, even exploring the U(l) and 
translation symmetries of the quantum chain we are restricted to lattice sizes up to L = 24, 
where the largest sector we can calculate the lowest eigenenergy has dimension 5,136,935 
(L = 24, n = 13). Since the calculations are done for a fixed density of particles p the larger 
number of sites we can use, and consequently get a better precision for our estimates, is at 
the density p = 1/2. 
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1.1771 


-1.4540 


1.001 


0.169 


0.167 


1.490 


0.168 
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-0.99597 


1.1300 


-1.4425 


1.001 


0.169 


0.167 


1.490 


0.168 


1.50 


-0.9812 


1.0596 


-1.4295 


1.01 


0.169 


0.167 


1.48 


0.168 



TABLE II: The quantities e^u^c, x+x,Oj x —i,o an< l x o.i f° r the Hamiltonian (TTJ) with density 
p = 1/2 and some values of t p . The first line (t p = 0.75) corresponds to a point inside region 2 
and the second one (t p = 1) to a point at the line separating regions 2 and 3. The remaining lines 
correspond to points in region 3 (see Fig. [I]). 

In table [U we present some of our estimated values of several quantities that characterizes 
the critical behaviour of the quantum chain ([1]) inside region 3 (t p = 1.15, 1.15, 1.35 and 1.5). 
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The calculations were done at density p = 1/2. In the first two lines, for comparison, we 
also give the results for a point inside region 2 (t p = 0.75) and a point at the line separating 
regions 2 and 3 (t p = 1). We include in the table the estimated values for the ground-state 
energy per particle in the bulk limit (we set h = in ([I])), the sound velocity v s and the 
magnetic field h that fixes the ground state at density p = 1/2. 

The value of the exponent x+i,o = x~i,o = x p for t p = 0.75 and t p = 1 are known from 
the solutions of (USD, flHD and ([26]), i. e., x p = 0.16770 (tp = 0.75) and x p = 0.1666 (t p = 1). 
The comparison of these last results with the first two lines of table [TT] indicates that the 
errors are in the last digit. Like regions 1 and 2, region 3 is also massless with a conformal 
central charge c = 1. The dimensions x + i j0 and x_i >0 in table [III are associated to the sectors 
with densities P = | + x and p = \ — respectively. These dimensions in a standard 
Coulomb gas phase, like in region 1 and 2, are equal and correspond to the dimensions x±i t o 
in (El]). We also show in table [TT] the dimension xo,i related to the first gap with zero (mod. 
7r) momentum in the sector containing the ground state. The dimensions x±i$, as shown in 
the table, is close to l/4xo,i, in agreement with ( )24l) . The dimensions in the first two lines 
of table II (t p = 0.75 and t p = 1) although close are not equal, as we can also check in Fig H] 
Table II also indicate that inside region 3 the conformal dimensions almost does not change 
as we change t p . Since our results are valid only for the density p = 1/2 we do not know if 
this small (or no) variation remains valid for other points inside region 3. 

Once in region 3 of Fig. [1] by increasing the value of t p , for a fixed density, we reach region 
4. In this region several eigenlevels crossings happen in the finite lattice. In Fig. [6] we show 
for L — 10 (left) and L = 16 (right) these crossings. The crossings among eigenlevels with 
distinct momenta are visualized directly in the figure. The kinks in the curves are due to 
the level crossings with the same momentum. 

As we cross from region 3 to region 4 the ground-state energy shows a discontinuity on its 
derivative due to a change of its relative position with an excited eigenstate. Unfortunately 
also in this region, due to numerical instabilities, it is quite difficult to solve directly the 
BAE ([8]). The direct diagonalization of the quantum chain, except at the density p = 1/2, 
can only be done for quite few lattice sizes. Our estimate of the line separating regions 3 
and 4, shown in Fig. [1] is then just schematic. This indicates a first order phase transition 
along the line separating these regions. 

Inside region 4 and for p = 1/2 we verify that the leading finite-size correction of the 
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FIG. 6: Lowest eigenenergy E(p) with momentum p, as a function of t p for the Hamiltonian (JT]) at 
density p = 1/2 and lattice sizes L = 10 (left) and L = 16 (right). We set /i = in the figure. The 
energies E(p) are the lowest eigenenergies with momentum p. The eigenlevels with p ^ 0, ir are 
doubled degenerated (momentum p and — p). The kinks on the curves are due to a level crossing 
among the two lowest eigenenergies with the same momentum. 

ground-state energy is not 0(1/ L 2 ) as expected in a massless conformaly invariant phase 
(see (TTTj) ) . This is illustrated in Fig. [7] where we show for some values of t p the ground-state 
energy as a function of 1/L 2 for the density p = 1/2. We see in this figure the distinct 
behaviour for the points belonging to regions 2 and 3 as compared with those of region 4 
(lower two curves). The oscillatory behaviour in Fig. [7] for the points inside region 4 is a 
consequence of the level crossings happening inside this region. As we can see in Fig. [6] tor 
t p ~ 1.65 many crossings happen with the crossing position dependent of the lattice size. 
As we change t p inside this region, for a given finite lattice L, the momentum of the ground 
state also changes. For example for t p w 2.25 (see fig. [6]) the ground state for the quantum 
chains with L = 10 and L = 16 has a non zero momentum and is doubled degenerated. 
The oscillatory behaviour shown in Fig. [7J turns out the finite-size scaling analysis quite 
imprecise. Although not conclusive these oscillatory behaviour in energy and momentum 
indicate that we have in region 4 effects of incommensurability of the density distributions, 



similar as happening in other models 19] 



In region 5, as in region 3, it is difficult to solve numerically the BAE due to the occurrence 
of complex roots mixed with real ones. In this case our analysis was restricted to small lattice 
sizes (L < 24). Since we should consider sequence of lattices with fixed density and on this 
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FIG. 7: Ground state energy per particle as a function of L~ 2 for the Hamiltonian ([T]) with h = 
at density p = 1/2. The lattice sizes are L = 6, 8, . . . , 24. The curves correspond to t p = 0.25, 0.5 
(region 1), t p = 0.75, 1.00 (region 2), t p = 1.25 (region 3) and tp = 1.85, 1.95 (region 4). 

region p > 1/2, the number of data we can have on a given finite-size sequence is quite small. 
However, for a given lattice size and density, as we change t p we did not see the crossings 
of eigenlevels observed in region 4 (see figure E]). This indicates that region 5, similarly as 
regions 1,2 and 3 is a critical c = 1 Coulomb gas phase. We greatly welcome more precise 
and convincing results for regions 4 and 5. 



IV. SUMMARY AND CONCLUSIONS 



We have made a detailed analysis of the spectral properties and phase diagram of the one 
of the new spin one model introduced in |6|]. This model is exact integrable for arbitrary 
values of the parameter t p . Our analysis was restricted to the cases where t p > 0. At t p = 1 
we have shown that the model recovers the deformed SU(3) Perk-Schultz model (or spin-1 
Sutherland model) at the special value of its deformation parameter q = e j27r//3 and external 
magnetic field. We verified that at this point the low lying eigenvalues of the model are 
the same as the corresponding ones of the XXZ quantum chain with anisotropy A = 1/2. 
We can then interpret the model we studied as an integrable generalization of the deformed 
SU(3) Perk-Schultz model or the XXZ quantum chain at anisotropy A = —(q + 1/q) with 
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q — gi27r/3 

The BAE, whose solutions give the eigenspectrum of the model are quite difficult to solve 
analytically or numerically for general values of the parameter t p and density of particles 
(magnetization). Our results are summarized in Fig. (TJ where we have regions 1-5. 

Regions 1,2,3 and 5 belong to a critical phase governed by an underlying Coulomb gas 
conformal field theory with critical exponents varying continuously. We distinguished these 
regions according to the BAE roots of the low lying eigenenergies of the quantum chain. 
In regions 1 and 2 these roots are real. This fact allowed us to solve directly the BAE for 
quite large lattices (L ~ 1000). Exploring conformal invariance we obtain good estimates for 
the conformal anomaly and anomalous dimensions of operators of the underlying conformal 
field theory. In these two regions we could take the thermodynamic limit and obtained the 
critical exponents in terms of integral equations (see (|25l)-(l26l)). 

In regions 3,4 and 5 the BAE roots corresponding to the ground state are not real and 
very difficult to calculate even for relatively small lattice sites. In these regions our analysis 
were based on the direct calculation of the eigenspectra for lattices sizes L < 24. 

Our results indicated that, regions 3 and 5, although having BAE complex roots in the 
ground state, have the same critical behaviour as in regions 1 and 2. Actually regions 1,2,3 
and 5 are quite similar to the XXZ quantum chain in the presence of a magnetic field, with 
anomalous dimensions xi t7n given by (1241) . with the value of x p depending on t p and p. 

As we cross from region 3 to region 4 there is a discontinuity of the ground-state energy. 
This is due to a crossing of the two lowest eigenenergies, and we expect that regions 3 and 
4 are separated by a first order phase transition. 

Inside region 4 we found an oscillatory behaviour for the energy gaps and momentum 
of the ground state, as we change the lattice size or the anisotropy parameter t p . These 
oscillations are due to a large number of crossing of eigenlevels. These crossings made our 
finite-size analysis imprecise. We believe that probably such oscillatory behaviour is due to 
incommensurability effects on the charge (local magnetization) distribution in the lattice, 



19] . 



as happens in other models whose incommensurability is established 

We conclude mentioning two interesting open problems for the future. The derivation of 
the i?-matrix associated to these new integrable spin one model and the extension of the 
present study to the second exact integrable model introduced in p. This second model at 
t p = 1 is related to the XXZ quantum chain with anisotropy A = —1/2. This last model has 
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remarkable properties. It is related to the problem of enumerating alternated sign matrices 



201 ] and the Hamiltonian with open boundaries is the evolution operator of a conformal 



invariant stochastic model, namely, the raise and peel model [2l|. We then expect that the 



second model in 



6[ with values of t p 7^ 1 may also have interesting connections with other 



interesting problems in physics and combinatorics. 
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